Density of states for dirty <i-wave superconductors: 
A unified and dual approach for different types of disorder 
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A two-parameter field theoretical representation is given of a 2-dimensional dirty d-wave supercon- 
ductor that interpolates between the Gaussian limit of uncorrelated weak disorder and the unitary 
limit of a dilute concentration of resonant scatterers. It is argued that a duality holds between 
these two regimes from which follows that a linearly vanishing density of states in the Gaussian 
limit transforms into a diverging one in the unitary limit arbitrarily close to the Fermi energy. 

PACS numbers: 72.15.Rn, 71.23.-k, 71.55.-I 
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Extrinsic defects in high T c materials have observable 
consequences for the local and bulk density of states 
(DOS) of quasi-particles. This is seen by scanning tun- 
nelling microscopy nuclear magnetic resonance 
and neutron scattering Q . Popular wisdom distinguishes 
"weak" extrinsic disorder introduced by defects above or 
below the copper-oxide planes through oxygen doping, 
from "strong" extrinsic disorder caused by substitutions 
of in-plane copper atoms with ions such as zinc. 

Given that the symmetry of the superconducting state 
is predominantly of the d-wave type in high T c materials, 
it is thus not surprising that a considerable theoretical ef- 
fort has been devoted to the effect of quenched disorder 
on a d-w&ve order parameter by reducing the problem to 
two models. The first one is a model of Anderson local- 
ization for a single "nodal" quasi-particle described by 
a 2-dimensional (2d) non-magnetic random Bogoliubov- 
de-Gennes (BdG) Hamiltonian [^[-^6|. The second one 
is a Kondo model in which a single magnetic moment 
couples to a bath of 2d "nodal" quasi-particles jnj. At 
this level of approximation no attempt is made to de- 
scribe either the existence of quasi-particle interactions 
or the self-consistent effect of defects on the supercon- 
ducting order parameter, although these effects might be 
important (p| . In spite of the simplicity of either models 
many important theoretical issues remain unresolved to 
this date. 

Even in the simpler case of Anderson localization, the 
effect of disorder on the quasi-particle DOS of a d-wave 
superconducting state is poorly understood. The physics 
of localization is expected to dominate at long distances 
and small energies in 2d ■ This fact together with 

the existence of particle-hole symmetry below T c makes 
it essential to account for the level repulsion upon ap- 
proaching the Fermi energy ep. However, the physics 
of level repulsion and localization are non-perturbative 
with respect to standard perturbative techniques [expan- 
sions in (fcp^)" 1 <C 1, fcp the Fermi momentum and £ the 
mean free path in the Born approximation] fljsfl . The 
existence of four isolated nodes of the superconducting 
pairing function Ak in the 2d Brillouin zone also inval- 



idates, upon approaching £p, expansions in powers of 
(k-pi)^ 1 <C 1 0. Moreover, there are some theoretical 
arguments and numerical evidences that have been at- 
tributed to the existence of nodes of Ak suggesting that 
the dependence on energy of the DOS depends in a dra- 
matic fashion on the microscopic details of the disorder. 
On the one hand, for a high density of very weak s-wave 
scatterers [the Gaussian limit (GL)] level repulsion dom- 
inates and causes the DOS to vanish algebraically upon 
approaching £p with an exponent that can depend on the 
disorder strength if the disorder is sufficiently smooth, 
i.e., correlated in space d]|,[l(|. O n the other hand, a 
dilute density of uncorrelated and resonant s-wave scat- 
terers [the unitary limit (UL)] is argued to result in a 
diverging DOS |llfll ^Jf"! . The purpose of this paper is 
to present a unified picture of these two limits. 

Starting from a BdG Hamiltonian describing a dirty d- 
wave superconductor with white-noise correlated disorder 
that depends on two parameters, the impurity concentra- 
tion p and the variance g of the disorder, we construct 
an effective field theory that generates disorder averaged 
Green functions. We show how this effective action re- 
duces to the nonlinear sigma model (NLSM) describing a 
random BdG Hamiltonian with time reversal symmetry 
and spin rotation symmetry (class CI in the classification 
of Ref. in the limit of p — > oo with g finite. The ef- 
fective action is treated by performing an expansion in 
powers of p (virial expansion) in the opposite limit of g 
fine-tuned to its resonant value g — > g r with p finite. Car- 
rying out the expansion up to first (second) order in p 
allows to identify this limit with the UL (the expansion 
parameter). As was the case with the (fcp^) -1 expansion, 
the virial expansion converges only for not too small en- 
ergies. In order to better understand whether there can 
or cannot be a crossover from the UL to the GL, we 
use a variational method and show that the form of the 
self-consistent equations is controlled by how the chemi- 
cal potential is introduced in the problem. This suggests 
that the energy alone cannot connect the two limits if the 
chemical potential is chosen at the "resonant" condition 
for the UL. By expanding perturbatively for large ener- 
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gies, we identify a duality relation between the GL and 
UL series that would imply, at low energies, that the di- 
verging DOS l/\e\ In 2 \l/e\ is connected by duality to the 
vanishing DOS |e| in the GL arbitrarily close to e-p. This 
non-perturbative argument relies in an essential way on 
the nodes of At- 

To model a dirty 2d d-wave superconductor, we use 
the tight-binding Hamiltonian Ti = Hq + V, 

n » = - E *« E c Usa + ^ 4Ai + Hc - ' w 
v = E E ^ c L^, 

that describes electrons hopping on a square lattice made 
of A/" sites with lattice spacing a. The real valued hopping 
matrix tjj contains in its diagonal the chemical poten- 
tial fiQ. Its off diagonal elements describe nearest neigh- 
bor hopping (t), next-nearest hopping (t') etc., and yield 
the normal state dispersion in reciprocal space. The 
order parameter Ay = (Ao/2) (6ij±jc — $i,j±y) enforces 
the d x 2_ y 2 symmetry. The disorder potential Vi describes 
N uncorrelated s-wave scatterers with the individual po- 
tential strength Vq. Moreover, N is randomly distributed 
according to a Poisson distribution. Hence, 
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where the overline denotes disorder averaging (the sub- 
script c refers to the cumulant). Instead of parametrizing 
the disorder by (ft, Vq) it is more convenient to introduce 
the parametrization (p, g) where g = pV 2 . We can then 
identify two distinct limits: 



p — > co, g finite. GL 
g — > g r , p finite. UL 



(2b) 
(2c) 



The resonant disorder strength g r = oo when the normal 
state dispersion is particle-hole symmetric, whereas it is 
a finite energy scale that depends on pq, t' ', etc., other- 
wise. To simplify notation we will assume a particle-hole 
symmetric normal state dispersion from now on: g T = oo. 

Effective field theories in Anderson localization are 
constructed to compute the disorder average over prod- 
ucts of the Green functions G(z) = (iz — Ti)^ 1 at differ- 
ent Matsubara energies z. For the problem at hand, the 
(bosonic) replicated canonical partition function 

— 1 i N — l 

V[<j>]e--^( i!1 - n °) ( P[A{p,g)f , (3) 



a=l , 



yields, upon insertion of a source term, the disorder av- 
erage of the single-particle Green function G(z) whereby 
disorder averaging is defined by weighing equally all the 
spatial configurations that N impurities can take on the 
lattice. Here, 7i = Ti-o + V is a 2Mn x 2JVn matrix 
representation of Eq. ([!]) (in the Nambu representation 
with the Pauli matrices 71,2,3 to represent the particle- 
hole grading) and is a replicated bosonic real field with 
2xJVxn components. In order to make a connection be- 
tween the GL and the UL we need to relax the condition 
that N is fixed in Eq. ([5j). This is done by defining 
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This Poisson grand canonical partition function is nor- 
malized to one in the replica limit n — ► 0. Denoting av- 
eraging with Z am by ((• • •)), we deduce that the average 
number of impurities is (N) = a and that the standard 
deviation in the number of particles scales like ^[u. From 
now on we will identify the impurity concentration p with 
a/M. The interacting field theory described by Eq. (^) 
has already been introduced by Mezard, Parisi and Zee 
as an effective field theory for random Euclidean matri- 
ces Jl9| (see also Refs. pp|-p2]|). Our derivation shows 
that the interacting field theory with a Gaussian interac- 
tion appears very naturally in the problem of Anderson 
localization with Poisson distributed impurities. 

According to Eq. (^) , the GL for a white noise distri- 
bution of impurities should be recovered as p — > 00. It 
is natural to expand the exponential of S\ in this limit. 
The first non-vanishing contribution to the action in the 
grand canonical partition function is of order g 



Zln = J V[<t>] exp(-Sf-Sf), 



!(V T 7 3 </>) , (5) 



aside from a shift in the chemical potential Sp = J pg. 
Equation (jsj) describes the effective action needed to cal- 
culate the average Green function of the random BdG 
Hamiltonian H.% + v. The superscript "g" indicates that 
the chemical potential /Iq = p a — ^fpg in Hq. The ran- 
dom potential v is white noise correlated and Gaussian 
distributed with vanishing mean and variance g/A, re- 
spectively. The usual justification for neglecting terms 
of order p(g/p) 3 ^ 2 is that they are irrelevant at the fixed 
point to which S^ + Sf flows. This fixed point is itself de- 
scribed by a NLSM on the manifold Sp(2n) in the limit 
of weak disorder, i.e., when the mean free path in the 



2 



Born approximation ~ 2t/g, is much larger than the lat- 
tice spacing a and for energies close to the band center 
f4 = pp5|. At this fixed point the DOS v{e) is only 



sensitive to the existence of nodes for not too small ener- 
gies |2^,Q • The physics of weak localization comes into 
play below the energy scale ef ~ (B/A) cxp(— it B 2 /2g) 
at which a pseudo gap opens up. Below an even smaller 
energy scale ef ~ ef [g/(A 2 B 2 )] exp(— A 2 ) random ma- 
trix theory (RMT) predicts that the DOS becomes lin- 
ear and insensitive to the nodes of Ak (the prediction of 
RMT can break down for correlated disorder in which 
case sensitivity of the DOS to the nodes of Ak is recov- 
ered ||jjie|). Here, B cx y/WA^ is the band width and 
A oc y^2i/Ao is the Dirac cone anisotropy. For notational 
simplicity 2t — Ao = 1 from now on. 

According to Eq. ( pcj ) the natural expansion for small 
impurity concentrations is the virial expansion , i.e., one 
in powers of p of cxp(— Si) in Eq. (Q). We have calcu- 
lated the increment (5z/ p )(e) in the DOS due to p = 1 
and p = 2 impurities, respectively. As it should be we 
reproduce the known increment 5v^ (e) induced by one 
impurity off and at resonance (g = g T ) to first order in 
p p^MRH . Off resonance the DOS crosses over to a 
linearly vanishing DOS upon approaching ep. This is 
not so at resonance. We infer from the calculation of 
the averaged single particle Green function (G(z)) up to 
second order in p that the virial expansion is an expan- 
sion in the small parameter a u = Trp/[z 2 ln(z 2 )]. Hence, 
the virial expansion is not uniform in the Matsubara en- 
ergy z. This non-uniformity of the virial expansion is 
the counterpart to the fact that in the GL the expansion 
in X/kp£ of (G(z)) is really an expansion in the small 
parameter a g = g\n(z 2 )/n 0). [The simplest guess, as 
suggested by the first terms in the virial expansion, for 
a small expansion parameter that interpolates between 
these two limits is a = a u a s /(a u + a s ).] The non- 
uniformity of the perturbative expansions in a s,u reflects 
the importance of quantum interferences already at the 
level of the single-particle Green function. Both pertur- 
bative breakdowns below the thresholds ef' u are signa- 
tures of the existence of isolated nodes of Ak in the 2d 
Brillouin zone that require non-perturbative approaches 
to access the DOS arbitrarily close to £p. We have not 
been able to find a counterpart to the Sp(2n) NLSM that 
describes the GL and is amenable to a calculation of the 
DOS arbitrarily close to £p. On the other hand, compar- 
ison of the expansions in a e and a n up to second order 
suggests that the GL and UL are related by g/ir <-» np 
and ln(z 2 ) <-*■ l/z 2 ln(z 2 ). If this relationship holds to 
all orders and extends to the non-perturbative regime 
as well it would allow access to lim e _> ep v(e) in the UL 
from knowledge of lim e ^ EF v{e) in the GL. We will argue 
from a detailed study of a variational approach that this 
duality relating the GL and UL indeed holds. 

For practical purposes ef' u can be considered small 
and the crossover of the DOS between the GL and UL 



is needed when e > maxef' u . This can be achieved 
by treating the interacting theory (^) within a varia- 
tional approximation. We thus replace the kernel of 
So + Si + pN in Eq. (||) by the non-interacting, transla- 
tionally invariant, and replica diagonal kernel iz — TL SC , 
z' = z + z. In reciprocal space ?i k f = £ k 7 3 + Aj c 7 1 , 
e k = e k — /io + £> and A k = A k + A. The self-consistent 
parameters z, e, and A, i.e., the disorder induced self- 
consistent shifts in the self-energy, chemical potential, 
and gap, respectively, are obtained by solving the self- 
consistent equations 



n(z') = p 

£{*') = p- 

z l 

A(z') = p 



z 



z 2 + e 2 ■ 
e 



A 2: 



A 2 



A 



z 2 + e 2 + A 2 ' 
1 

A7 



(n,£,A)(,') = l£- 



At) 



k (z') 2 + (4) 2 + (A k ) 2 ' 



(6a) 
(6b) 
(6c) 
(6d) 



The self-consistent DOS is ^ sc (e) = lim^^o ^Hz^e+in 
Im [i sgn(^) n(z')] . It is sufficient to consider Eqs. ( |6a| , |6b| ) 
alone since A fex , fey = -A feyjfca; allows A = 0. 

That details of disorder matter in 2d <i-wave supercon- 
ductors (see Ref. jH| ) already follows from the variational 
approach. Indeed, Eqs. ( |6a| , |6b| ) are not so much sensitive 
to the energy z as they are to the chemical potential 
when it comes to describing the crossover from the UL to 
the GL. Equations ([ja 6b) reduce to the two well-known 
self-consistent equations 



p — > oo, g < oo, po 



Spg^n(z') 



z 

j 
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Po, 
(7a) 



p < oo, ^o = 



^Q(z') = ^, e = 0, 

z 



(7b) 



in both the GL and the UL HO- U is thus 

imperative to scale po with ^fpg for e to be tuned away 
from the UL resonance |e| -C \z\. 

The self-consistent Eq. (|7^) does not capture the inter- 
play between the physics of localization and that of level 
repulsion below the energy scale ef p,H- The same 
should be true of Eq. ((7^) below the energy scale e" at 
which the virial expansion breaks down. We are going 
to argue that the physics of localization and level repul- 
sion results in a diverging DOS upon approaching ep in 
the UL. Our argument relies first on identifying a dual- 
ity relating the perturbative expansions of the self-energy 
in the GL and UL valid at large energies relative to the 
disorder strength. We then postulate that the same du- 
ality extends to small energies relative to the disorder 
strength. Let z be a large Matsubara energy, i.e., g (p) 
is sufficiently small so that level repulsion can be ignored 
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and the physics of localization is that of the standard or- 
thogonal symmetry class. In this regime, we can safely 
remove the self-consistency condition in 



n(z + z) = J2\ £ 

p=o I P 



4 Mz>)] +P (sgn)(z), gl, 



C [n(z')]- p (8Sil)(z), UL, 



(8a) 



by substituting fl(z + gQ(z + ...)•••) [Q,(z + p/Q(z + 
...)•••)] with Q{z) + (d z Q)(z)n(z) g + 0(g 2 ) [Q(z) + 
(d z lad)(z) p + 0(p 2 )]. In this way, the self-energy re- 
duces to the perturbative expansion 

o ( , + g) .iME!) + ^ 2 ) + 2](^ f 

77 n { 2z> ln(^) 

_2zy, r/3|[ln(^ 2 )] +p (f) P , GL. 

Here, the combinatoric coefHcient /3|' u depends on the 
order in g, p to which the self-consistency has been trun- 
cated. For example, the substitution Sl(z') = Q(z) gives 
Pp = 1/p(p — 1) = /9p an d accounts only for the pure rain- 
bow diagrams in the expansion of the Dyson equation for 
the self-energy in the GL. Truncation to higher order in g 
accounts for diagrams that are not purely of the rainbow 
type but in which impurity lines never cross. Neglecting 
diagrams with crossing of impurity lines is only justified 
for sufficiently weak disorder at a fixed energy Q . Cor- 
respondingly, the condition on the disorder strength for 
the physics of localization and level repulsion to be unim- 
portant is 



-a g = ^-ln(l/z 2 ) < 1 

7T p 



Z > Z, 



- a 



z 2 \n(l/z 2 ) 



< 1 



oh 



z > z, 



GL. 



Ob 



UL. 



(9a) 
(9b) 



Note that the expansions in the GL and UL are related 

(10) 



9 , / ix 1 

-^tt P , \n(z ) <-> 2TT2V 
7T z z m{z' : ) 



By undressing the propagator in Eq. (jSq), two expansions 
of the same scaling form Q(z')—Cl(z) — 2(z/n){^[\n(z' 2 )+ 
2]a s ' u — /(a g,u )} have been derived for large energies 
in the GL and UL, respectively. We now assume that 
the self-energy can also be written as some scaling func- 
tion / that can be expanded in powers of a s ' u outside of 
the radius of convergences (^), i.e., at low energies. The 
scaling function / should be specified by the symmetry 
class CI, i.e., the physics of level repulsion induced by 
the particle-hole symmetry of a time-reversal and spin- 
rotation symmetric superconductor. This scaling func- 
tion is not completely unknown since RMT predicts that 
/(a s ) tx a s for |a g | > 1. [(-2z/7r)( 5 lnz 2 /7r) is the 
combination that, upon analytical continuation, leads to 
v s {e) oc |e| arbitrarily close to e — 0.] By the duality 



assumption f(a u ) oc a u must also hold outside of the 
radius of convergences ([)!]) [leading to the combination 
(— 2z/-k){ 2 T p 3 )]. Upon analytical continuation, dual- 
ity then predicts a diverging DOS v n {e) oc l/|e|ln |e| 
arbitrarily close to the Fermi energy £ F in the UL. The 
robustness of this conjecture requires a careful study of 
the role played by diagrams with crossing impurity lines 
(the Goldstone modes in NLSM terminology) in the UL 
that is beyond the scope of this paper. 

Our main conclusion is that the crossover between the 
resonant UL and the GL in a 2d dirty d-wave supercon- 
ductor is not controlled by energy only but by the chem- 
ical potential and the strength of the disorder, which to- 
gether tune out of the resonance. A variational approach 
yields self-consistent equations that interpolate between 
the GL and UL. In practice these self-consistent equa- 
tions describe accurately the DOS since the energy scales 
signalling breakdown of perturbative expansions are ex- 
pected to be rather small. We have established under 
what conditions the diverging DOS l/\s\ In 2 |e| predicted 
by Pepin and Lee in Ref. can arise. A duality be- 
tween the UL and the GL that holds in the regime of 
applicability of the variational approach must extend to 
the non-perturbative regime. In this way the apparent 
dichotomy between a DOS controlled solely by RMT in 
the GL and one controlled solely by the nodes of the gap 
is elegantly resolved. Only the leading functional depen- 
dence on energy is predicted to be universal by RMT. 
Accordingly, we cannot predict the prefactor of the di- 
verging DOS whereas Pepin and Lee state that this pref- 
actor is exactly given by p, coincidentaly the prefactor 
for independent impurities. 
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